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Within nuclear thermodynamics and an isospin-dependent transport model we 
investigate respective roles of the nuclear mean field and the 2-body stochastic 
scattering on the evolution of density and isospin fluctuations in either mechan- 
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Q\ ■ ically or chemically unstable regions of neutron-rich matter. It is found that the 
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^ ■ the 2-body scattering influences significantly the later growth of the isospin fluc- 



mean field dominates overwhelmingly the fast growth of both fluctuations, while 



"q ■ tuation only. Moreover, both fluctuations grow in mechanically unstable systems, 

while only the density fluctuation grows significantly in chemically unstable ones. 



Furthermore, the magnitude of both fluctuations decreases with the increasing 



isospin asymmetry because of the larger reduction of the attractive isoscalar 
mean field by the stronger repuslive neutron symmetry potential in the more 
neutron-rich matter. Finally, several experimental measurements are proposed 
to test these findings. 
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The rapid advance in experiments with rare isotopes has opened up several new frontiers 
in nuclear science The planned Rare Isotope Accelerator (RIA)will further enhance the 
exploration of these forefronts dramatically 0. Besides detailed information on the structure 
of exotic nuclei in many unexplored regions of the periodic chart, it is now possible to study 
novel properties of isospin-asymmetric nuclear matter with extreme neutron to proton ratios. 
Prospects for discovering new physics in warm neutron-rich matter of various densities that 
can be created transiently in nuclear reactions induced by exotic neutron-rich nuclei have 
generated much interest in the nuclear science community In particular, the density 

dependence of the symmetry term in the equation of state (EOS) has recently received 
much attention since it is among the most important but very poorly known properties of 
neutron-rich matter [p~0| pT5[] . This term is very important to the mechanisms of Type II 



supernova explosions and neutron-star mergers. It also determines the proton fraction and 
electron chemical potential in neutron stars at f3 equilibrium. These quantities consequently 
influence the cooling rate of protoneutron stars and the possibility of kaon condensation in 
dense stellar matter [|T^-|2D|. 



Based on the EOS obtained from various microscopic many-body theories and phe- 
nomenological models PTJ-OBJ, it has long been predicted that isospin-asymmetric nuclear 



matter under certain conditions can be mechanically or chemically unstable, i.e., 
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where P, p n and 5 = (p n — p p ) / (p n + p p ) are the pressure, neutron chemical potential and 
isospin asymmetry, respectively. In the mechanically or chemically unstable regions, small 
fluctuations in density and/or isopin- asymmetry are generally expected to grow. However, 
information about the respective growth rates of these fluctuations, how and why they 
may depend on the isospin-asymmtry of nuclear matter, seeds of fluctuations and their 



isospin-dependence, how and which aspects of nuclear dynamics cause the growth of these 
fluctuations is rather rare. This information is vitally important for understanding the 
structure and stability of both neuton stars and radioactive nuclei as well as mechanisms 
of nuclear multifragmentation in reactions induced by neutron-rich nuclei. In this Rapid 
Communication, we obtain such information within nuclear thermodynamics and an isospin- 
dependent transport model. We also propose several experimental measurements to test our 
findings. 

The boundaries of the mechanical and chemical instabilities in the configuration space 
of density p, isospin asymmetry 5 and temperature T can be established readily by using a 



thermal model for isospin asymmetric nuclear matter |24| . We use here a Skyrme-type phe- 
nomenological EOS for isospin-asymmetric nuclear matter in the parabolic approximation 



27,28]. The energy per nucleon at T = can be written as 



e(p, 5) = °- + -^—u° + \eW" + So(po) • « 7 • S 2 . (3) 
2 l + O" 5 

In the above u = p/po is the reduced density, e% is the Fermi energy in symmetric nuclear 
matter at normal density, and a = —123.6 MeV, b = 70 .4 MeV and a = 2 corresponding 
to a stiff nuclear EOS of isospin-symmetric nuclear matter. The last term is the symmetry 
energy whose density dependence is currently rather uncertain pl|-|r5|j. We adopt here a 
parameterization used by Heiselberg and Hjorth- Jensen in their recent studies of neutron 
stars with S (p ) = 30 MeV and 7 as a free parameter. We note that a value of about 
7 = 0.6 was obtained by fitting to the result of variational many-body calculations PP| , |3T)| . 
The nucleon chemical potential in asymmetric nuclear matter p q ,(q=n for neutrons, q=p for 



protons) at temperature T is thus |3T 



p q = au + bu a + v ^ sy + T 



ln(^i) + V Vl±l hn{ ^ 



(4) 



■> 1/ 2 

where = 2wh 2 / (m q T) is the thermal wavelength of a nucleon. The coefficients b n are 



obtained from mathematical inversion of the Fermi distribution function ||31|| . The single 
particle symmetry potential v ^ sy is 
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where "+" and "-" are for neutrons and protons, respectively. The corresponding pressure 
is 



P = lap u 2 + b^-u° +1 + (So 7 Po« 7+1 

L (7 + 1 



Z n=l 

We then use the Gibbs-Duhem relation 
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to find boundaries of mechanically unstable regions (ITS: isothermal spinodal) in the p — 5 
plane for given temperatures. Since the chemical instability condition has to be evaluated 
at constant pressures, the following Maxwellian relation 
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T,P \ / T,p \ / T,8 \ / T,5 \ / T,p 

is used to find boundaries of the chemically unstable regions (DS: diffusive spinodal). 

Shown in Fig. 1 are the boundaries of the mechanical (thick lines) and chemical (thin 
lines) instabilities in the p — 5 plane with 7 = 0.5 at T = 5 (lower window), 10 (middle 
window) and 15 MeV (upper window), respectively. It is seen that the diffusive spinodal 
extends further out into the plane and envelopes the region of mechanical instability; the two 
regions of instability do no overlap. As the temperature increases, both instabilities become 
less prominent over a more narrow range of densities at smaller isospin asymmetries. These 
features are independent of the parameters of the EOS and are in good agreement with 
those based on more microscopic many-body theories |f2~T|-f2"3|1 . In particular, we found that 
the variation of 7 parameter has very little effect on the instability boundaries. Therefore, 
in the following a constant of 7 = 0.5 is used. The results shown in Fig. 1 not only provide 
the motivations but also serve as a guidance for our following numerical simulations. 

The evolution of isospin-asymmetric nuclear matter in the mechanically or chemically 
unstable regions is studied by using the isospin-dependent transport model . The model 



uses consistently the isospin-dependent EOS and the corresponding potentials as in the 
thermal model outlined above. Moreover, scatterings among neutrons and protons are fully 
isospin-dependent in terms of their total and differential cross sections as well as the Pauli 
blockings. Nucleons are initialized in momentum space using a Boltzmann distribution 
function in a cubic box of length L box with periodic boundary conditions. We have checked 
that the use of Fermi-Dirac distribution does not alter our results at all as we are mainly 
interested in the evolution of reduced fluctuations with respect to the initial state which 
is unstable against numerical fluctuations in the mechanical or chemical instability region. 
Thus, the initialization with the Boltzmann distribution is sufficient and much more efficient 
numerically. The box is further divided into cells of 1 fm 3 volume in which the average 
density p ce u and isopsin asymmetry S ce u are evaluated. We used 10 4 test particles per 
nucleon in evaluating the p ce u and 5 ce u. Shown in Fig. 2 is an illustration of the evolution 
of a system initialized in the mechanically unstable region with = 5 MeV, Si = 0.6 and 
Pi = 0.05 fm~ 3 . The most interesting feature shown here is the gradually increasing isospin 
fractionation. This is indicated by the spreading of the initial system into regions with p < Pi 
and 5 > 5i and where p > p« but 5 < Si. The variations of p and S with respect to their initial 
values can be characterized quantitatively by using, respectively, cr d (t) = (p 2 — p 2 ) 1 ^ 2 and 
<Js(t) = (5 2 — S 2 ) 1 ^ 2 , where the average is over all cells. Furthermore, the degree of isospin 
fractionation can be quantified by using the ratio (N/Z) gas /(N/Z)n quid , where (N/Z) gas and 
(N/Z)n quid is the isospin asymmetry of the low (p/po < 1/8) and high (p/po > 1/8) density 
region, respectively. 

Shown in Fig. 3 are the reduced variation with respect to the initial state in 
isospin asymmetry as(t) — cr^ (0) , density <Jd{t)/Pi and the degree of isospin fractionation 
(N j Z) gas /(N/Z) u qu id as a function of time for a system initialized at Tj = 5 MeV, p, =0.05 
fmT 3 and Si =0.2, 0.6 and 0.9, respectively. As a reference and check of our approach, re- 
sults (dash-dot lines) are also shown for a system initialized with Tj = 15 MeV and Si = 0.9 
where it is both mechanically and chemically stable. For this system, it is seen that both the 
isospin and density fluctuations stay almost constant and there is no isospin fractionation at 
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all. While for the system initialized in the mechanically (chemically) unstable region with 
5i =0.2 and 0.6 (Si =0.9), it is seen that both the isospin and density fluctuations grow 
faster with the decreasing isospin asymmetry Si, and the opposite trend is observed for the 
strength of isospin fractionation. Moreover, both the density and isospin fluctuations grow 
in the mechanically unstable systems, while only the density fluctuation grows significantly 
in the chemically unstable ones. It is also interesting to note that the isospin fractionation 
happens later in the more neutron-rich matter. In the latter early fluctuations are smaller 
and thus take longer time to grow before they can trigger the isospin fractionation. 

We now proceed to explore the seeds of fluctuations and their isospin dependence, and 
to investigate how and which aspects of nuclear dynamics are important in governing the 
growth of fluctuations. We shall also study the question why the more neutron-rich matter 
is more stable against both the density and isospin fluctuations. The evolution dynamics 
of asymmetric nuclear matter is governed by the isospin-dependent nuclear mean filed and 
stochastic nuclear scatterings. In our approach, first-order effects of 2-body stochastic scat- 
terings are included through the collision integral of the BUU equation. Therefore, there are 
two main seeds of fluctuations in our approach, i.e., the early numerical fluctuations from 
the random sampling of the initial state and the later 2-body stochastic nuclear scatterings. 
Both of them may lead to the growth of fluctuations by propagating through the nuclear 
mean field in the early and later stages of the evolution, respectively. Shown in the upper 
window of Fig. 4 are the average number of collisions per nucleon as a function of time for 
S = 0.2 and 0.9, respectively. Before about 40 fm/c, there is essentially no collision because 
of the strong Pauli blocking when the phase space nonuniformity due to the initial fluctuation 
is still small. Later, nuclear scatterings become important, moreover, they are more frequent 
with the decreasing isospin asymmetry. The less frequent 2-body scatterings observed with 
the higher S is because of the isospin-dependence of both the nucleon-nucleon cross sections 
and the Pauli blocking rates. It is well known that the cross section for neutron-neutron 
scatterings is only about 1/3 of that for neutron-proton collisions at beam energies below 
about 1 GeV f32[. One also expects from Fermi statistics that neutron-neutron scatterings 

6 



are more strongly Pauli blocked in the more neutron-rich matter. To investigate the respec- 
tive roles of the nuclear mean field and the 2-body scatterings, model studies by turnning off 
the nucleon-nucleon collisions have been performed. As shown in the middle window of Fig. 
4, the collisional seeds of fluctuations lead to the significant growth of the isospin fluctuation 
in the later stage of the evolution. However, they have very little effect on the growth of 
density fluctuations as shown in the lower window. A comparison of the results obtained 
with and without the collision integral indicates that the fluctuations are overwhelmingly 
dominated by the nuclear mean field. The observed isospin-dependence of the fluctuations 
can be understood from the interplay between the attractive isoscalar mean field and the 
repulsive symmetry potential for neutrons. The symmetry potential, as shown in Eq. [|, is 
repulsive for neutrons and attractive for protons and their magnitudes increase with the 
increasing isospin asymmetry. Thus, the resultant attractive mean field is weaker in the 
more neutron-rich matter. Because of the small number of scatterings, particularly in the 
early stage of the evolution, the growth of fluctuations is thus mainly determined by the 
strength and sign of the resultant nuclear mean field according to the linear response theory 
for asymmetry nuclear matter PS|,R3,R3. Based on the latter, the magnitude of both fluc- 



tuations can grow larger with the increasing strength of the attractive resultant mean field, 
and thus also with the decreasing isospin asymmetry 5. 

What are the important physical implications of our findings? We expect several effects 
that can be tested experimentally. It is well known that more neutron-rich systems are less 
bound and have smaller saturation densities. However, fluctuations and their growth are also 
important for determining the final state of a neutron-rich system, such as in the projectile 



fragmentation in producing exotic beams [35]. Our results above indicate that fluctuations 



actually have a compensating role to the lower binding energy in stablizing neutron-rich sys- 
tem where fluctuations are smaller and do not grow as fast as in symmetric ones. Moreover, 
our results on the isospin fractionation accompanying the evolution of fluctuations indi- 
cate that the configuration of a more dense, isopin-symmetric region surrounded by a more 
isospin-asymmetric gas as in halo nuclei is a natural result of the isospin-dependent nuclear 



dynamics. Furthermore, we expect that the multifragmentation and isospin fractionation in 
nuclear reactions induced by neutron-rich nuclei to happen on longer time scales compared 
to symmetric reactions of the same masses. These expectations can be tested by measuring 
products of multifragmentation, in particular, the neutron-neutron, proton-proton as well 
as fragment-fragment correlation functions, in comparative studies of isospin symmetric and 



asymmetric nuclear reactions |35|-p8 



In summary, we have investigated the isospin dependence of mechanical and chemical in- 
stabilities in neutron-rich matter within nuclear thermodynamics and the isospin-dependent 
transport model. We found that the mean field dominates overwhelmingly the fast growth 
of both density and isospin fluctuations, while the 2-body scatterings influence significantly 
the later, slower growth of the isospin fluctuation only. Moreover, both of the fluctuations 
grow in the mechanically unstable systems, while only the density fluctuation grows signif- 
icantly in the chemically unstable ones. Furthermore, the magnitude of both fluctuations 
decreases with the increasing isospin asymmetry because of the larger reduction of the at- 
tractive isoscalar nuclear mean field by the stronger repuslive neutron symmetry potential 
in the more neutron-rich matter. Finally, several experimental measurements are proposed 
to test these findings. 
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FIG. 1. Boundaries of the mechanical (thick lines) and chemical (thin lines) in the den- 



sity-isospin asymmetry plane at a temperature of 5, 10 and 15 MeV, respectively. 
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FIG. 2. An illustration of the evolution of asymmetric nuclear matter in the density-isospin 
asymmetry plane. Each point in the scatter plots represent one cell of 1 fm 3 volume in a cubic 
box of side length Li, ox = 10 fm. 
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FIG. 3. Evolution of the reduced isospin fluctuation (upper window), reduced density fluctu- 
ation (middle window) and strength of isospin fractionation (lower window)in a cubic box of side 
length Lb ox = 30 fm and density 0.05 /m -3 . The dash-dot lines are calculated with Tj = 15 MeV 
and 5i = 0.9; while the solid, dot and dashed lines are calculated with Tj = 5 MeV and 6i =0.2, 
0.6 and 0.9, respectively. 
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FIG. 4. The average number of successful nucleon-nucleon collisions per nucleon (upper win- 
dow), the isospin (middle window) and density (lower window) fluctuations as a function of time. 
The initial conditions are the same as in Fig. 3. The dot and dash-dot lines are results of calcula- 
tions without the 2-body scatterings. For comparisions on the same scale, a factor of 7 is multiplied 
to the isospin fluctuations with 5 = 0.9 in the middle window. 
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